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Abstract. - We investigate the non-Markovian quantum dynamics of a hybrid open system con- 
sisting of one qubit and one qutrit by employing a stochastic Schrodinger equation to generate 
diffusive quantum trajectories. We have established an exact quantum state diffusion (QSD) equa- 
tion for the dissipative qubit-qutrit system coupled to a bosonic heat bath at zero temperature. 
As an important application, the non-Markovian QSD equation is employed to simulate the en- 
tanglement decay and generation measured by negativity. Finally, some steady state properties 
of the hybrid system are also discussed. 
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Introduction. — The recent development of quantum 
information and quantum optics [UH] has triggered much 
effort devoted to physics of qubit systems in various phys- 
ical settings as they are the most basic building blocks 
of quantum computing and quantum information process- 
ing. The dynamical behaviors of open quantum systems 
(qubit, qutrit, etc) are abundantly documented for the 
cases where the environmental noise is weak and the band- 
width of the noise is broad such that the so-called Born- 
Markov approximation gives rise to a reliable description 
' of the dynamics of the quantum open system [3] . Under 
this approximation, one may employ either the standard 
Lindblad master equations or their alternatives includ- 
ing quantum trajectories, Fokker-Planck equations and 
Langevin equations, to name a few [3H3- When the ex- 
ternal environment is a structured medium or the system- 
environment coupling is strong, then the non-Markovian 
equations of the motion such as non-Markovian master 
equations have to be used to describe the time evolution 
of the reduced density operator of the system of interest 

[SHE]. 

A stochastic Schrodinger equation named non- 
Markovian quantum state diffusion equation by Strunz 
and his coworkers has been applied to many interesting 
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physical systems including a dissipative two-level sys- 
tem, quantum Brownian motion model and a multiple 
qubit system [T3H23] . A recent application of the non- 
Markovian QSD to molecular aggregates can be found 
in [53]. It has been known that the higher dimensional 
quantum systems (qutrit, qudit, etc.) have demonstrated 
potential applications in quantum information processing 
25, 26 . Their dynamical aspect in a non-Markovian 
regime is still not well understood due to the lack of exact 
non-Markovian QSD equation or exact master equation. 
In the framework of quantum open system, a qubit-qutrit 
system is clearly the first non-trivial extension beyond 
the qubit-qubit system, yet it is sophisticated enough to 
provide useful insight into non-Markovian dynamics of 
high dimensional open systems [TUHT]- I n the case of 
entanglement dynamics, the qubit-qutrit system allows 
a rigorous characterization of entanglement when the 
negativity is employed. It should be noted that the 
disentanglement pathways for a qubit-qutrit system have 
been studied when the qubit and qutrit are coupled to 
individual Markov environments, respectively [28], but a 
systematic investigation into the qubit-qubtrit model in a 
non-Markovian regime is still missing. 

Purpose of this Letter is to generalize the non- 
Markovian QSD approach to the case of the qubit-qutrit 
system. In particular, we have derived the exact time-local 
quantum state diffusion equation for the qubit-qutrit sys- 
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tern, in which the explicit formation of O-operator is very 
different from that in the two-qubit case [35]. In what 
follows, we will first introduce the general formalism of 
non-Markovian quantum trajectory. We then derive the 
exact time-local QSD equation for the qubit-qutrit dissipa- 
tive model at zero temperature. Finally, by employing the 
non-Markovian QSD equation, we investigate the negativ- 
ity dynamics of the hybrid system for different anisotropic 
coupling parameters, environmental memories and initial 
states. The details for the exact QSD equation are left to 
the two appendices. 

General formalism. — The model considered in this 
paper is described by the qubit-qutrit system coupled lin- 
early to a general bosonic environment consisting of a set 
of harmonic oscillators: 
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where H sys = uj a S^ + uj b S b and L = S_ + kS b . k 
depicts the asymmetry between the coupling coefficients 
of the two subsystems with the environment, and explicitly 
the notations are given by 
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In what follows we will show that a time-local stochas- 
tic Schrodinger equation at zero temperature called QSD 
equation can be derived from the above microscopic 
model. Using the Bargmann state bases for the environ- 
ment degrees of freedom, the total quantum state includ- 
ing the system and the environment could be written as 
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. = \\z 1 )\\z 2 )---\\z x ) with \\z x ) = e**°i|0), 
z\ 2 = Ea M 2 , and [&] = ... ^ .... It was 

shown that the wave function of the central system i/jt(z) 
is governed by a linear stochastic differential equation with 
a functional derivative of the pure state of system over the 
stochastic noise z [T41IT51: 
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-iH sys + Lz t — L' I dsa(t, s) 
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where a(t, s) is the environmental correlation function and 
z t — —iY^\9\ z \ e%wxt i s a complex Gaussian process sat- 
isfying M[z t ] = M[z t z s ] = and M[z*z s ] = a(t,s). The 
density matrix of the system can be recovered by the en- 
semble average of quantum trajectories: 

p t =M[\Mz))(Mz)\}= I —e-^ 2 \Mz))(Mz)l (5) 



One of the advantages of the non-Markovian quantum tra- 
jectory method is that a small number of trajectories could 
also be used to estimate qualitatively the dynamical be- 
haviors of the open system [37] . Crucial to application of 
the general QSD equation ([J} to a practical problem is to 
transform eq. (@| into a time-local stochastic differential 
equation: ^j^- = 0(t, s, z)ij: t {z)- Then we have 
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Mz) = [-iH sys + Lz t - L+0(t, z)]Mz), (6) 



where 0(t, z) is defined as J Q a(t, s)0(t, s, z)ds. Although 
there is no general recipe to construct the O-operator ex- 
plicitly, the exact O-operator for this qubit-qutrit dissipa- 
tive model can be constructed by the following expression: 
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1, 2, • • • , 13, are given by 
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where a 5. are the Pauli matrices. And the coefficient 



functions /j's, j — 1,2, ana Pj's, j 

could be determined by the consistency condition |13j : 
s _ aw* it turns out to be 
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As a primary result in this paper, we have explicitly de- 
rived the O-operator for the qubit-qutrit model, therefore 
we have established for the first time the exact dynamical 
equation for this hybrid system. From the initial condition 
for the O-operator 0(s,s,z) — L, we have fi(s, s) = 0; 
fj(s,s) = 1, j = 2,3,4,5; fj(s,s) = k, j = 6,7,8; 
Pj(s, s, si) = 0, j = 9, 10, 11, 12; and pi 3 (s, s, si, s 2 ) = 0. 
The partial differential equations and the boundary con- 
ditions for fj's and p^s are given in Appendix A. 

The linear QSD equation in eq. ^ does not conserve 
the norm of the state vector ipt(z), so for numerical simu- 
lations, one typically employs the nonlinear version of the 
QSD equation © [15]: 



j t ^ t {z) = [-iH sys + A t (L)z + A t (tfO(t,z)) 
- (L^) t A t (0(t,S))}Mz), 



(9) 
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where the norm-preserved wave function is defined as 
Mz) = Note that z t = z t + f Q a*{t,s){tf) s ds 

is the shift noise, A t (Q) = Q — (Q)t for any operator Q, 
and (Q)t = (V'tlQIV't) denotes the quantum average. Ex- 
plicitly, the operator 0(t, z) may be written as 

8 12 . t 

0{t,£) = J2 F i& D 3 + 12 P j (t,8i)z^d8iD j 
3 = 1 j=9 Jo 

+ Pi3(t,si,s 2 )z Sl z S2 dsids 2 Dx 3 , (10) 

Jo Jo 

where the definition of Fj(t)'s, j = l,---,8, Pj(t,Si)'s, 
j = 9, • • • , 12, and Pi3(i, Si, S2) can be found in Appendix 
A. 






K=0 




■— -K=0.25 




---k=0.5 


K=1.0 











5 

cot 



10 



Fig. 1: The negativity as function of cjt with (a) different 7's 
(k = 1), and (b) different «'s (7 = 0.3) obtained by an av- 
erage over 1000 trajectories. The initial state is chosen as 
(1/V2)(|1)a|1)s + |0)a|0)s). 
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Fig. 2: The negativity as function of cjt with (a) different 7's 
(k = 1), and (b) different k's (7 — 0.3) obtained by an average 
over 1000 trajectories. The initial state is chosen as |2}a|0)s. 

Simulations and Discussions. — Entanglement dy- 
namics of the qubit-qubit system in Markov and non- 
Markovian baths has been extensively studied in the last 
decade [30M40], Equipped with the exact QSD equation 
for the qubit-qutrit system, we are in the position to study 
rigorously the non-Markovian entanglement and decoher- 
ence dynamics for this hybrid system. For the purpose 



of entanglement dynamics, qubit-qubit (2 ® 2) and qubit- 
qutrit (2 <g> 3) systems are particularly interesting due to 
the fact that entanglement measured negativity [4TH44] 
serves both sufficient and necessary conditions for an ar- 
bitrary state of these two bipartite systems. Precisely, the 
negativity is defined as N(p) = \\p Tn \\ — 1 = 2^ IMj'Ij 
where Tb stands for the partial transpose with respect to 
one of subsystems B, and fij's are the negative eigenval- 
ues of the matrix p Tu . A qubit-qutrit state represented 
by the density matrix p is entangled if and only if N(p) is 
positive. The negativity N(p) varies from zero for all sep- 
arable states to unity for the maximally entangled states. 
For a general pure state of qubit-qubit and qubit-qutrit 
system, it reduces to the definition of Wootters' concur- 
rence and generalized concurrence for pure states [45 . 46 . 
respectively. 

Once the correlation function a(t, s) is chosen, it is 
straightforward to calculate these coefficient functions in 
eq. (fLOj) . In the following, for the sake of simplicity, we con- 
sider a Gaussian complex noise z t satisfying the Ornstein- 
Uhlenbeck process: a(t,s) = ^e~ 7 '*~ s ', where T is the 
dissipation rate and 7 describes the memory time of the 
environment. More precisely, y -1 defines the finite cor- 
rection time of the environment. When 7 — > 00, the cor- 
relation function approaches the Markov limit: a(t,s) = 
TS(t — s). In Appendix B, the details of the differential 
equations for all the functions Fj's (j = 1, • • • , 8) and P/s 
(j = 9, • ■ ■ , 13) are provided. It should be noted that the 
Pj terms contain explicitly the history integral over com- 
plex Gaussian noise z t . It is expected that they will be 
important in non-Markovian dynamics when 7 is not too 
large (far from Markov regimes). 

In the following numerical simulations, we choose loa = 
ujb = ui, T = 1 and the double integral term about Pi 3 
containing noises is neglected to save the computer mem- 
ory. It is easy to see that the higher orders are clearly less 
important in a weakly non-Markovian regime. 

In figs, [lja) and []Jb), the initial state is a maximal 
entangled state |^(0)) = (1/V2)(|1)a|1)b + |0)a|0) b ). 
Fig. [lja) shows the effect of 7 on the decay rates of 
the entanglement of the system. When 7 is as large as 
3, the dynamics is seen to be very close to the Markov 
limit. It shows the reservoir can quickly ruin the coher- 
ence of the initial state, which evolves into a separable 
state monotonously. When 7 = 0.1, the environment has 
a long memory time, so that the non-Markovian features 
become dominant, hence the negativity exhibits a strong 
fluctuation pattern and maintains a highly entangled state 
for a long time. When 7 = 0.3, it shows a moderate 
non-Markovian behavior. Fig. [ljb) illustrates how the 
anisotropic coupling parameter k affects the entanglement 
evolution. The asymmetrical coupling slows down the en- 
tanglement decay gradually as the parameter k becomes 
smaller and smaller. The extreme case is that k — 0, 
where one of the subsystems is completely decoupled from 
the influence of the environment. 

Figure. [U[a) plots the entanglement generation of the 
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separable state |^(0)) = 12)^10)5. If the qubit is not 
interacting with the qutrit, then the entanglement gen- 
eration is purely due to the indirect interaction between 
the subsystems induced by the common bath and memory 
times. For this initially separable state, the larger 7 (near- 
Markov reservoir) leads to a faster transition from |2) to 
|1) compared to the smaller 7 (a stronger non-Markovian 
regime). Interestingly, the near-Markov reservoir in turn 
causes a quicker generation of entanglement with one ex- 
citation transition in a short time interval. For some ini- 
tially separable states, the non-Markovian property may 
be essential in entanglement generation. In the case of 
7 = 0.1, representing a highly non-Markovian regime, it 
clearly shows that the dynamics generates a higher degree 
of entanglement than that in the case of 7 = 0.3 after 
tut w 9. While fig. ^h) describes the effect of different 
k's on the entanglement generation with the same 7. It is 
very interesting to see that the entangling power is highly 
sensitive to the coupling balance for the qubit and qutrit 
interacting with the common bath. 
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Fig. 3: The distribution of the ground state p\\ (black solid 
lines) and the negativity of the steady state N s (blue dashed 
lines) as functions of k with different initial states in the long- 
time limit. 

It is interesting to examine the long-time behavior of the 
combined qubit-qutrit system under the influence of the 
non-Markovian environment. For the Ornstein-Uhlenbeck 
noise, the information for the long-time limit may be ob- 
tained by solving the Markov master equation for the 
qubit-qutrit system 

d tPt = -i[H sys ,p t ] + ^[L,p t tf] + j[Lp u tf]. (11) 

Now we consider several interesting initial states in the 
four subspaces spanned by (i) |2)^|1)b, (ii) 12)^10)^ and 
\l) A \l) B , (in) |1)a|0) b and |0)a|1)b, (iv) |0)a|0) Bi respec- 
tively. When T = 1, oja = wb, we have 



dtPee 



-(1 + n 2 )p 66 , 



0tP55 = K 2 p 66 - - (p 45 + p 54 ) 



P55, 



dtP54 



K P66 - 7:(P44 + P55) ~ (1 + Tr)P54, 



9tP33 = P55 - P33 - -^(P23 + P32) + n(p i5 + p 5i ) + K 2 p 4i , 

9tP32 = P5i - + y)P32 ~ -j(P22 + P33) + K/?44, 
dtP22 = P44 - 2 (P32 + P23) - K 2 p 2 2, 



dtPll 



P33 + n(P23 + P32) + K 2 p 2 2, 



where the notations for the six basis vectors are |6) = 
\2) a \1)b, |5) = \2) A \0) B , |4) = ||3) = \1)a\0)b, 

|2) = |0)a|1)b, and |1) = \0) a \0}b- From the above equa- 
tions, one can easily conclude that for the stationary so- 
lutions we always have pee = P55 = P54 = P44 = (they 
are independent on the initial states). But P33 = — K/932, 
P22 = -P32/K) Pn = 1 - P33 - P22, which are dependent 
on the initial states. Therefore, we have 



pu(oo)\00) A B\ 
[l-pn(oo)]|V(K))<V(«)|, 



(12) 



dtpa = Pee - (1 + K 2 )p 4i - -(/O54 + P45), 



where |V(«)> - (1/VT+^){k\1)a\0)b - |0>a|1)b), which 
is an eigenstate of the Lindblad operator L = + kS 1 ^. 

The final value of pn is a function of the k and the 
initial states as well as the negativity of the final state. 
The black-solid line with the square marker in Fig. ([3]) 
describes the probability distribution pn of the ground 
state |0)a|0) b for the initial state |V>(0)) = \2} A \0} B - In 
the range < K ^ 1, it decreases with k monotonously. 
And the negativity [equals to 2(l — pn)n/ (1 + k 2 )] is shown 
by the blue-dashed line with the same marker. It coincides 
with the curves in Fig. [2jb). We also plot the numerical 
results (the two lines with the circle marker) for the initial 
state \</>(k)) = (1/VT+~k?)(k\2)a\0)b - \1)a\1)b), which 
is also an eigenstate of L, but this one is not protected 
by the common dissipative environment. It shows clearly 
that the state evolves into a mixture of the ground state 
|0)a|0)b and \i/j(k)) given by (fT2)) . An interesting result 
arising from the above discussions is that the residue en- 
tanglement degree reaches its maximal value at K ~ 0.75 
instead of 1 which corresponds to the balanced coupling. 
It is also interesting to see that the final state could be- 
come an approximate pure state when k rs 0.55, which 
demonstrates a rather different patterns from those ini- 
tially separable states considered in this section. 

Conclusion. — In summery, we have investigated 
non-Markovian QSD approach to the qubit-qutrit sys- 
tem coupled to a zero-temperature bath. With the non- 
Markovian QSD equation, the entanglement dynamics of 
the qubit-qutrit system is studied carefully. The numer- 
ical simulations have been used to show the negativity 
dynamics for different initial states. We also discuss the 
long-time behaviors of the initially entangled or separable 
states of the hybrid system. Entanglement generation and 
decoherence have been shown to be sensitively dependent 
on coupling strength, memory time and initial states. The 
approach described here provides a fundamental tool for 
qubit-qutrit systems when the non-Markovian properties 
become important. 
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Appendix A: Equations of the Motions for the boundary conditions are: 

Coefficients and Boundary Conditions - General 

Case. — By the time evolution equation for O-operator 



p 9 (t,s,t) = -Kfi(t,s) + f 2 (t,s) - f 3 (t,s), 



in eq. (|8]), we get a set of partial differential equations Pw(t, s, t) — nfi(t, s) + f±(t, s) — f$(t, s), 

about the coefficient functions in eq. 0: Pll (t, s, t) = nf 2 (t, s) - nf 4 (t, s) + f 6 (t, s) - f 7 (t, s), 



dtfi(t,s) = i{2u A -u B )fi + F 1 f 3 + F A fi + KF A f 3 
- nFsfi - KP w (t, s), (Al) 



pi 2 (t,s,t) = Kf 3 (t,s) - nf 5 (t,s) + h{t,s) - f s (t,s), 
p 13 (t,s,t,si) = 0.5[np 9 (t, s,si) - Kp w (t, s,si) 

+ p n (t,s,sx) -p 12 (t,s,sx)]. (A14) 



dtf 2 (t s) — iloa} 2 + F 2 f 2 — Fife — F 3 f 2 — KFi/e Appendix B: Equations of the Motion for the 

F F P (t P ( (A2) ^ oe ^ clen ^ s: Ornstein-Uhlenbeck Noise. — In the 

+k 6/2 k 7/2 9 ( ,s) k u( ,s), ( ) case f the Ornstein-Uhlenbeck process, by the definitions 

and the boundary conditions in eq. (|A14|) . we get a set of 
d t f 3 (t, s) = iuj A h + F 3 f 3 - nF 5 f 7 + F 7 fx + nF 7 f 3 closed differential equations for F/s (j = 1, ■ • • , 8), P,'s 
- KF s f3-nP 12 (t,s), (A3) (7 = 9,-,12)andJ^: 

dtF^t) = (-^ + 2iw A -iw B )F 1 +F 1 F3+F 1 F A 
d t h{t,s) = ioj A U + F 1 f 7 + F4/4 + F4/7 - F 5 / 4 - KFiFa + KFaFi-KFiFs-KPio, (Bl) 



F s h-P 10 (t,s), (A4) 



d t F 2 (t) = -^ + {-j + i ( JA)F 2 + F 2 2 -F 1 F 6 -F 2 F 3 



d t f 5 (t,s)=iu>Af 5 + F 5 f 5 + K F 5 f 8 , (A5) ' ' 2 

+ kF 2 F 6 - kF 2 F 7 - K FiF 6 - P 9 - kP u ,(B2) 

dtf6(t,s) = iuj B f& + F 2 f 6 - F 4 f 6 + nF 6 f 6 

- F 7 f 2 -P n (t,s), (A6) 9 t F 3 (t) = - 1 l + (~~ f + iL J A)F3 + F 1 F 7 + F 3 2 + K F 3 F 7 

- kF 3 F 8 - kF 5 F 7 - kP 12 , (B3) 



d t f 7 (t,s) = iLo B f 7 + F 3 f 7 -F 5 f 7 + F 7 U + K F 7 f 7 
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- F 8 f 3 - P 12 (t,s), (A7) d t F 4 (t) = -^- + {- 1 + iLu A )F i + F 1 F 7 -F 1 F % +Fl 

d t f 8 (t,s)=iu> B f 8 +F s f 5 + KF s f 8 , (A8) - F 4 F 5 + F 4 F 7 -P 10 , (B4) 

a u , . ^ P „ , „ d t F 5 (t) = ^ + (- 1 + iu A )F 5 +Fg + KF 5 F 8 , (B5) 

d t p 9 (t, s, si) = 2ioj A p 9 + F 2 p 9 - nF 5 pn + nF 6 p 9 I 



-KF 8 p 9 + P9/3 - nP 10 f 6 + P 1 1/1 



+ K P 11 / 3 - K P 1 2/2-2 K P 13 (^ S , Sl ), (A9) Wi) = — + (-7 + ^)P 6 + ^6 

- F 4 P 6 + kF 6 2 -Pii, (B6) 

d t pio(t,s,si) = 2iu> A Pw + Ftfv + F^pxo + nF4p 12 

+ P0/5 + kPio/s, (A10) 5 t P 7 (i) - ^ + (- 7 + J wb)F7 + P3P7-P 3 P8 

+ F4P7 - P5P7 + «P 7 2 - P12, (B7) 

+ K F 6 p n -F 8 p 9+ P 9 f 7 -P w f e + Pnh d t F s (t) = ^ + (-j + ^ B )F S + F 5 F 8 + K Fi , (B8) 
+ KP 11 f 7 -P 12 f 2 -2P 13 (t,s,s 1 ), (All) 

a t P 9 - (-2 7 + 2^)P 9 + ^(-kPi+P2-P3) + PiPii 

d t px 2 (t, s, si) = i(wa +w s )pi2 +P3P12 +P7P10 + F 2 P 9 - kF 2 P 12 + P3P9 + KP3P11 - KF5P11 

+ KF 7Pl2 + P 12 f 5 + K P 12 f s , (A12) + kP 6 P 9 - kP 6 P 10 - kF 8 P 9 - 2kP 13 , (B9) 

p-5 



Jim Jingl,2^ Ting Yul^ 



d t P 10 = (-2 1 + 2ilo a )P 10 + ^( K F 1 +F 4 -F 5 )+F 1 P 12 1 15 
+ F 4 P 10 + kF 4 P X2 + F 5 P 10 + kF 8 P 10 , 



(BIO) [16 

r . I 17 

dtPll = (-27 + %W A + iUB)PLl + ^ {kF 2 - kF± + F & t 18 

- F 7 ) + F 2 P n -F 2 P 12 + F i P n -F 5 P 11 -F 6 P 10 [19 

+ K F 6 P 11 +F 7 P 9 + K F 7 P 11 -F 8 P 9 -2P 13 ,(B11) [ 20 

[21 

22 

<9 t Pi 2 = (~2-f + iujA + iuB)Pi2 + ^-{KF 3 -KF 5 + F 7 

- F 8 )+F 3 P 12 + F 5 P 12 + F 7 P 10 [24 



«F 7 Pi 2 +/tF 8 Pi 2) 



(B12) 



d t Pi 



T7, 



[25 
[26 



't-ri3 = + 2zUM + iuj B )Pl3 + -J-(K P 9 ~ kF W 

+ Pu - P12) + F 2 P 13 + F a P 13 + kF 6 Aa B 
+ P9P2 + /sF 8 As + A0A1 + kPuPw, (B13) [29 



where Pj = P,(i) = J* a(t,s)P l (t,s)ds, j = 9, • • • , 12; [30] 

P13 = /o Jo s i) a (A s )Pi3(A s j si)dsdsi; And the ini- 
tial conditions for Fj(0),j = 1, • • • , 8, P/(0), j = 9, • • • , 12 
and Pi3(0) are all zero. 
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